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ABSTRACT 

If  a large  negative  voltage  is  applied  to  an  insulated  disc  inserted 
in  a collision-dominated  plasma,  the  electrons  are  expelled  from  the  region 
near  the  disc  leaving  a bare  cloud  of  ions.  This  "transient  sheath"  is  well 
documented  in  ion  acoustic  wave  experiments.  This  paper  examines  the 
evolution  of  the  ions  in  this  sheath  into  a collision-dominated  plasma.  A 
new  sheath  whose  dimension  is  of  the  order  of  twice  the  transient  sheath 
dimension  is  found. 
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EXPANSION  OF  AN  ION  SHEATH  INTO  A COLLISION -DOMINATED  PLASMA 

K.  E.  Lonngren^ 

The  study  of  sheaths  surrounding  biased  electrodes  inserted  in  plasmas  is  a well  studied 
problem.  For  steady  state  plasmas,  we  can  easily  cite  the  work  of  Langmuir  as  being  representa- 
tive. A sheath  which  was  later  discovered  by  Alexeff  and  his  coworkers  was  the  transient  plasma 
sheath*  *.  In  this  sheath,  a large  negative  potential  was  assumed  to  be  instantaneously  applied 
to  an  electrode  (say  a large  flat  plate)  in  the  plasma.  Owing  to  the  mass  difference,  the  electrons 
are  blown  into  the  plasma  leaving  a bare  cloud  of  ions  in  the  region  0 < x < x^  . The  plate  with 
voltage  -<J>0  applied  to  it  is  located  at  x = 0.  At  the  transient  sheath  dimension  x^,  the  potential 
and  the  electric  field  are  assumed  to  be  zero.  Solving  Poisson's  equation  in  the  transient  sheath 
region  with  the  above  boundary  conditions,  the  dimension  of  the  sheath  xT  can  be  computed  to  be 
Kj  = [2  e$0/ne2/t0]J  in  MKS  units. 

As  time  progresses,  the  ions  in  the  transient  sheath  are  "eaten"  by  the  electrode  and  eventu- 
ally a steady  state  Langmuir  sheath  is  set  up.  In  the  Langmuir  sheath,  there  exists  a balance 
between  the  amblpolar  ion  current  to  the  electrode  and  the  space-charge  limited  current.  For  time 
scales  between  the  creation  of  the  transient  sheath  and  the  eventual  formation  of  the  Langmuir 
sheath,  the  description  of  the  physics  involves  a numerical  solution  of  the  nonlinear  partial  dif- 
ferential equations  of  continuity  and  motion  coupled  with  Poisson's  equation  and  a Boltzmann 
relation  for  elections.  Widner,  et  al  have  done  such  a calculation  and  obtained  good  agreement 
with  experiment*2* 

In  that  study.  It  was  found  that  the  ions  within  the  transient  sheath  decayed  in  time  and 
one  could  approximate  the  ion  behavior  at  x = 0 as 


P+(t) 


1 

1+t 


(1) 


where  a suitable  normalization  has  been  Included.  Secondly  it  was  found  that  the  front  of  the 


sheath  expanded  into  the  plasma,  initially  at  a velocity  faster  than  the  ion  acoustic  velocity  c . 
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At  a distance  equal  to  approximately  2 xT  from  the  electrode,  it  was  found  that  the  sheath 
expansion  slowed  down  to  the  ion  acoustic  velocity.  Third,  large  Ion  currents  were  initially 
drawn  by  the  electrode  as  the  ions  started  to  move.  A similar  result  was  reported  also  by 
others* Finally,  Oskam  and  his  colleagues  noted  from  experiment  that  the  time  that  it  took 
for  a sheath  to  reach  an  equilibrium  state  (say  a Langmuir  sheath)  decreased  if  the  plasma  density 
Increased  or  if  ions  of  a smaller  mass  were  used*^’7\ 

The  work  that  will  be  described  here  suggests  that  the  above  observations  can  be  under- 
stood within  the  framework  of  a model  which  assumes  that  the  ion  velocity  is  proportional  to  the 
electric  field,  i.e.  a mobility  model.  In  Section  II,  we  present  the  model  and  describe  the  sheath 
evolution  in  Section  III.  Section  IV  is  the  conclusion  which  relates  the  predicted  behavior  with 

l 

these  numerical  and  laboratory  experimental  results. 

II.  The  Model. 

We  shall  assume  tha*  .he  plasma  can  be  described  by  the  set  of  one-dimensional  equa- 
tions 


8t  = 

8x 

(2) 

8E  _ 

(3) 

*0  8x 

P+  - P. 

1 = 

8x 

p+  iT  = w>+E 

(4) 

VX  B ~ 

CL+ 

II 

o 

(5) 

which  are  the  equation  of  continuity,  Poisson's  equation,  definition  of  current  and  Maxwell's 

equation  respectively.  The  set  of  equations  (2-5)  has  been  studied  by  Many  and  Rakavy  with 

particular  emphasis  on  the  theory  of  space-charge-limlted  currents  in  solids.'  ' Self  similar 

(9) 

solutions  for  the  set  have  recently  been  obtained.  De  Oliveira  and  Leal  Ferreira  have  studied 
this  set  with  reference  to  charge  neutralization  between  capacitor  plates  if  charge  of  one  sign 
filled  one  half  the  distance  between  the  plates  and  charge  of  the  other  sign  filled  the  other 
half.**0)  We  shall  make  extensive  use  of  the  philosophy  of  this  calculation  in  the  derivation 
that  follows. 
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We  find  the  characteristics  of  the  system  which  define  the  flow  lines  x(t)  for  the  positive 
mobile  charges 


= PE  (x(t),  t)  . 

The  acceleration  along  a flow  line  is  given  by 
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where  r = , p is  the  unperturbed  density  and  we've  used  (2),  (3),  (5)  and  (6).  Along  the 

*0 

flow  line,  the  variation  of  positive  charge  is 
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where  we've  used  (2), (4)  and  (6).  Equations  (7)  and  (8)  are  the  starting  equations  for  the  deriva- 
tion of  the  ion  sheath  evolution. 

Ill  Sheath  Evolution 

At  r = 0,  we  shall  assume  that  a large  negative  potential  -4>Q  is  applied  to  the  elec- 
trode at  x = 0.  Before  the  ions  can  move,  the  electrons  are  "blown"  Into  the  semi-infinite  uni- 
form plasma  to  at  least  a distance  of  the  transient  sheath  dimension  x,j,  . ^ The  temporal  be- 
havior of  the  ions  in  the  transient  sheath  region  0<  x < x^,  is  given  by  the  solution  of  (8)  with 
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As  time  Increases,  ions  will  be  "eaten"  by  the  electrode  and  will  expand  into  the  plasma. 


As  we  are  more  interested  in  the  behavior  for  x > x^  , we  shall  assume  an  insulating  coating 


covers  the  electrode  such  that  It  draws  no  current.  In  fact,  experiments  by  Chen  and  Schott 
confirm  that  this  technique  is  a very  efficient  method  to  launch  Ion  acoustic  waves. 

At  x:  x^,  the  electric  field  can  be  computed  from  Gauss  Law  to  be 

*T 

E(x_,t)  = PJT>  T-  (10) 

1 + *0 

The  velocity  of  ions  for  x > x^,  can  be  computed  from  (7)  to  be 
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The  time  rQ  Is  the  time  that  a particular  ion  leaves  x = x^  . The  density  ratio  p /pQ  is  the 

density  ratio  of  excess  electrons  that  are  "blown"  out  of  the  transient  sheath  into  the  neutral 

x p0 

plasma.  We  shall  assume  that  it  Is  uniform  for  1 < — < (1  + — ).  The  value 
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where  (10)  has  been  used.  We  note  from  (11)  that  the  velocity  of  the  Ion  cloud  slows  down 


time  increases. 


The  position  of  the  ions  in  the  plasma  can  be  determined  by  Integrating  (11)  from  rQ  for 


ions  at  the  sheath  edge  x^  to  x at  a time  r . We  find 
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The  leading  edge  of  the  ion  expansion  occurs  for  the  ions  that  leave  the  transient  sheath  at  tq  = 0. 

A new  sheath  criterion  is  determined  as  follows.  Let  us  find  the  time  T and  position  L 
where  the  leading  edge  of  the  ion  sheath  expansion  slows  down  to  the  ion  acoustic  velocity 
Cg  - Cf  « q/Vpq-  From  (9),  (11)  and  (12),thls  time  T is  defined  as 

P. 
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Combining  (13)  and  (14),  we  obtain 
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This  can  be  simplified  by  using  the  definitions  for  mobility,  the  Ion  acoustic:  velocity  and  the 
transient  sheath  to  yield 


= 1 
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where  w . Is  the  Ion  plasma  frequency  with  density  p and  v is  the  collision  frequency.  We 
pi  o 

note  that  If  the  excess  electrons  resided  In  the  region  x^  < x < 2x^  and  pQ  = p _,  then  L/x^  < 2. 

The  Ion  density  for  x > x^,  could  In  principal  be  found  by  eliminating  the  parameter  tq  be- 
tween (6),  (11),  (13)  and  Poisson's  equation.  As  this  cannot  be  carried  out  explicitly,  we  follow  the 

(10) 


procedure  given  by  de  Oliveira  and  Leal  Ferreira 

PoVp. 


Differentlate  (13)  with  respect  to 
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Along  the  flow  lines,  the  electric  field  can  be  found  from  (6),  (11)  and  (12)  to  be 

P 

(t-t) 

P0  0 

E = E(xt,  tq)  « 

Therefore  (13)  can  be  rewritten  using  (10)  and  (18)  as 

p</p-  *0  e(xt,to) 

XT 
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(18, 


(19) 


Now  consider  rQ  to  be  a function  of  x and  r and  differentiate  (19)  with  respect  to  x and 


use  (19)  and  Poisson's  equations.  One  finds  the  ion  density  p+(x  > x^,)  to  be 


P+<*  > 


/pn 


■^,t't«T  p-  1 


(20) 
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To  obtain  p+(x  > as  a function  of  x and  r,  we  assign  values  to  tq  and  use  equations  (13) 
and  (20).  Typical  results  for  the  value  p-/*0  =1  as  computed  with  an  HP-25  programmable 


hand  calculator  are  shown  In  Figure  1.  We  note  the  rapid  neutralization  of  the  charge 
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in  the  region  x > xT  . 

IV  Conclusion 

The  evolution  of  an  ion  sheath  into  a collision  dominated  plasma  has  been  computed.  The 
model  predicts  the  temporal  ion  behavior  in  the  transient  sheath  and  the  slow  down  of  the  sheath 
evolution  observed  by  Wldner,  et  al.  ^ The  unnormalized  time  that  the  sheath  reaches  an  equi- 

- P.  *P0  - P e - 

librium  (say  dx/dr  = C ) is  computed  from  (14)  as  — t = constant  = t 

3 . m p0  « EoV 

Therefore  this  equilibration  time  t « agrees  with  the  Oskam,  et  al  experimental  findings.  ' ' 
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Figure  1 


Figure  1. 


Normalized  ion  density  as  a function  of  distance  x/x^  and 

time.  Electron  density  in  plasma  >Jchosen  such  that 

p /pQ  = 1.  Note  the  slowdown  of  the  front  and  the  rapid  charge 

x 

neutralization  in  the  region : 1 < — — < 2 . 
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